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ABSTRACT 

Under  single  scattering  conditions,  water  droplets 
clouds  do  not  depolarize  the  backscattered  light. 
However,  backscattered  light  from  multiple 
scattering  will  be  depolarized.  The  level  of 
depolarization  is  a  function  of  the  droplets  size, 
the  cloud  extinction  coefficient  value  and  profile; 
it  has  also  an  important  dependency  on  the  lidar 
field-of-view  (FOV).  The  use  of  depolarization 
information  to  retrieve  cloud  microphysical 
properties,  using  Multiple-FOV  has  been  the 
object  of  studies,  [1],  [2].  Recently  the  use  of  the 
depolarization,  at  a  single  FOV,  has  been  studied 
for  cloud  with  linear  liquid  water  content  profiles, 
[3],  [4].  In  this  paper  we  present  the  mechanism 
leading  to  depolarization  and  identify  the  FOV 
values  for  which  the  information  on  particle  size  is 
high.  Also  Monte  Carlo  simulations  for  cloud 
with  constant  and  ramp  up  profiles  are  presented. 
The  degree  of  linear  depolarization  as  a  function 
of  cloud  penetration  is  significantly  different  for 
both  cloud  profiles.  This  suggests  that  the  use  of 
the  degree  of  linear  depolarization  at  a  single  FOV 
should  be  used  with  caution  to  determine  clouds 
micro-physical  parameters. 

1.  DEPOLARIZATION  RATIOS 

Backscattered  light  from  spherical  particles  does 
not  change  the  polarization  orientation  state  of  a 
linearly  polarized  incident  light.  This  particularity 
allows  the  discrimination  of  spherical  particles 
from  non-spherical  particles.  However,  under 
multiple  scattering  conditions,  even  spherical 
particles  generate  depolarized  signals.  The  linear 
depolarization  ratio  (5//n)  at  a  scattering  angle  ft  is 
defined  as  the  ratio  of  the  power  integrated  over 
the  azimuthal  angles  cp  in  the  plane  of  polarization 
orthogonal  (f±)  and  parallel  (I//)  to  the  linearly 
polarized  source: 


Using  Mie  theory,  it  can  be  shown  [1]  that: 
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Where  in  Pi,  P2  and  P3  are  the  Mie  scattering 
phase  matrix  elements  defined  by  Deirmendjian. 

Figures  1  shows  the  phase  functions  and  the 
depolarization  ratios  at  a  wavelength  of  0.532  pm 
for  linear  incident  waves  on  water-droplets  of 
diameters  4  pm  and  20  pm.  The  water  droplets 
size  distribution  is  log-normal  with  a  geometric 
standard  deviation  of  0.1;  the  distribution  was 
truncated  for  values  smaller  than  0.02  the 
maximum  of  the  distribution. 
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Figure  1.  Phase  functions  and  linear  depolarization 
ratios  for  4  and  20  pm  in  diameter  water  droplets. 

For  depolarization,  the  major  observations  are: 


•  The  light  is  not  depolarized  in  the  forward 
scattering  direction  for  angles  as  large  as 
20°. 

•  The  depolarization  is  equal  to  0  at  180°. 

•  The  depolarization  increases  very  sharply 
around  180°. 

Figure  2  shows  the  depolarization  as  a  function  of 
scattering  angles  ranging  from  170°  to  180°  for 
water  droplets  of  4,  8,  12,  16  and  20  pm  (same 
distribution  as  for  figure  1).  For  smaller  particles, 
the  depolarization  peak  shifts  toward  smaller 
scattering  angles.  However  the  maximum  values 
keep  very  close  to  each  other. 


Figure  2.  Depolarization  ratios  for  water  droplets  of  4, 
8,  12,  16  and  20  pm  in  diameter.  The  rectangular  box 
represents  a  1°  width. 

2.  SCATTERING  AND  DEPOLARIZATION 

Double  scattering  events  and  higher  scattering 
events  leading  to  depolarization  are  illustrated  in 
Figure  3.  For  clarity,  only  a  single  backscattering 
has  been  illustrated  at  distance  zc  from  the  lidar 
system.  This  backscattering  is  preceded  by  a 
forward  scattering  event  that  can  come  from 
different  distances  za  .  A  lidar  with  a  given  FOV  ( 
)  will  measure  the  backscattered  light  without 
any  distinction  regarding  the  forward  scattering 
event  position  in  space.  The  second  order 
scattering  signal  detected  will  be  the  sum  of  all 
those  signals  and  it  will  be  given  by: 
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where  P(zc )  =  P0  exp[— 2  J a(z)dz  ] 
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where,  C  5  T  and  A  are  respectively  the  speed 
of  light,  the  pulse  width  and  the  collecting  optics 
area;  ocs  is  the  scattering  extinction  coefficient; 
p(r,z  ,p )  an(j  p(r,zc,j3b  )  are  the  values  of 
the  phase  function  for  the  forward  ( P )  and 
backward  scattering  angles  ( =  n  -  /3  +  0  )  for  a 
particle  of  radius  r;  za  is  the  distance  of  the 
beginning  of  the  cloud;  zc  is  the  distance  where 
the  scattered  radiation  is  measured;  the  quantity 
fas(z)p(r,z  ,/3)l  represents  the  forward 
scattering  coefficient  while 

[as( zc  )p(r,zc,fib )]  represents  the 

backscattering  coefficient.  The  FOV  6  is  easily 
related  to  the  scattering  angle  P  in  the  simple 

R  _  zctanQ 

geometric  relation:  tan  P  -  - —  .  For 

compactness  of  the  equations,  the  particle  size 
dependence  of  the  functions  F  and  3lin  has  been 
dropped  out. 

The  depolarization  ratio  for  scattering  angle  within  J3t 
is  calculated  as  follows: 
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Figure  3.  Double  scattering  processes  and  higher 
scattering  processes  leading  to  backscattering 
depolarization. 


Figure  1  clearly  shows  that  the  depolarization 
ratios  for  scattering  angles  smaller  than  20°  are 
very  close  to  zero.  Considering  the  small  field  of 
view  covered  by  a  typical  lidar  detector,  the 
measured  depolarization  effect  must  come  from 
the  backscattering.  Figure  4  shows  the  function 
as  a  function  of  the  scattering 

angle  P  . 
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Figure  4.  The  product  of  the  forward  scattering  phase 
function  at  scattering  angle  (3  and  the  depolarization 
ratio  at  angle  f3'=  [3  -  n  as  a  function  of  the  scattering 
angle  [3. 

Clearly  the  product  of  the  forward  scattering  phase 
function  at  scattering  angle  (3  and  the 
depolarization  ratio  at  angle  (3  -  (3  -  71  as  a  function 
of  scattering  angle  (3  contains  information  on  the 
size  of  the  particle  and  a  Multiple  FOV  lidar  will 
be  able  to  extract  that  information  under  some 


conditions.  Can  a  single  FOV  depolarization 
measurement  provide  information  on  droplet  size? 
Is  there  an  optimal  FOV  to  obtain  information  on 
droplet  size? 

In  view  of  a  better  understanding  of  the  problem, 
Figure  5  shows  the  product  of  normalized 
cumulative  forward  scattering  and  depolarization 
ratio  as  a  function  of  scattering  angle.  The 
cumulative  signal  at  a  specific  scattering  angle  is 
the  signal  which  would  be  recorded  by  a  lidar  with 
that  specific  angle  as  FOV.  The  different  arrows 
are  for  different  cloud  geometries  and  lidar  FOVs. 
For  small  scattering  angles,  the  average  scattering 
angle,  <|3>  over  a  penetration  depth  Az  (Az  -zc  - 
0  n2zr+Az 

za )  is  equal  to:  <  p  >=  u — - - . 

Az 

From  Fig.  5  and  Table  1,  it  is  clear  that  for  average 
scattering  angles  greater  than  2.5  degrees  the 
product  of  cumulative  forward  scattering  and 
depolarization  ratio  reach  values  very  close  to 
each  other  (except  for  the  4  pm  droplets).  A  small 
FOV  would  thus  show  better  differentiation 
capability,  which  reinforces  the  interest  of  the 
depolarization  single  FOV  lidar  based  method,  but 
at  the  cost  of  an  often  significantly  reduced  signal. 


Figure  5.  Normalized  product  of  cumulative  forward 
scattering  and  depolarization  ratio  as  a  function  of 
scattering  angle.  The  different  identified  arrows  are  for 
different  cloud  geometries  and  lidar  FOVs. 

Table  1  Average  scattering  angle  within  lidar  FOV  of 
0.5  and  2  mrad  for  different  cloud  base  and 
cloud  penetration  depth. 
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FOV:  0.5 
(mrad) 

FOV:  2.0 
(mrad) 

Qoudbase(m)  Az(m) 

<-r>  (Deg) 

^R>(Deg) 

500 

10 

0.72  (Al) 

2.89 

1000 

10 

1.44  (A2) 

5.76 

2000 

10 

2.87  (A3) 

11.49 

500 

20 

0.37  (Bl) 

1.46 

1000 

20 

0.72  (B2) 

2.89 

2000 

20 

1.44  (B3) 

5.76 

500 

30 

0.25(a) 

0.98 

1000 

30 

0.48(Q) 

1.94 

2000 

30 

0.96(0) 

3.85 

So  far  we  have  focused  our  effort  on 
understanding  the  physics  of  the  depolarization 
effect  caused  by  multiple  scattering  and  its  relation 
to  the  lidar  FOV  and  did  not  address  the  effect  of 
the  cloud  extinction  profile.  In  the  following  and 
last  section,  the  effect  of  cloud  extinction  profile 
will  be  studied  using  Monte  Carlo  simulation. 

3.  MONTE  CARLO  SIMULATION 

The  Undique  Monte  Carlo  (MC)  simulator  is  a 
multithread  software  based  on  the  Bohren  and 
Huffman  Mie  scattering  computation  routine.  The 
simulator  is  made  to  reproduce  the  characteristics 
of  an  actual  Flash  Lidar  system.  It  consists  of  an 
emitter/receiver  system,  a  target  and  a  propagation 
range  including  aerosols  of  various  properties. 
The  particularity  of  the  Undique  MC  is  its 
capability  to  image  the  scattered  light  on  a  detector 
array.  See  [6]  for  details. 
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Figure  6  Degree  of  linear  polarization  as  a  function  of 
cloud  penetration  depth. 

The  different  parameters  for  the  MC  simulation 
are  listed;  Emitter:  532nm,  0.5  mrad  beam 
divergence,  Pulse  width  and  shape:  10  ns, 
Gaussian;  Receiver:  Receiver  diameter  and  focal 
length;  0.2m  and  0.527m;  Cloud  base:  1000m, 
depth:  50  m,  optical  depth:  2.5,  profiles:  Square 
and  Ramp  up. 

Figure  6  shows  the  degree  of  linear  polarization  as 
a  function  of  cloud  penetration.  Degree  of  linear 
polarization  was  used  to  facilitate  comparison  with 
[3].  The  cloud  extinction  profile  has  a  strong 
influence  on  the  degree  of  linear  polarization  (or 
depolarization  ratio)  and  such  an  effect  is  only 
partly  accounted  for  in  [3]  and  [4].  In  absence  of 
knowledge  of  cloud  extinction  profile  or  working 
hypothesis  [4],  the  simple  measurement  of 
depolarization  for  a  given  FOV  will  most  likely 
not  be  enough  to  provide  information  on  droplet 
size. 

CONCLUSION 

Droplet  size  information  is  contained  in  the  lidar 
depolarization  signal,  especially  at  small  FOVs 
where  the  higher  order  scattering  effects  are  less 
important  and  where  the  product  of  cumulative 
forward  scattering  and  depolarization  ratio  shows 
distinct  values  as  a  function  of  droplet  size  for 
constant  cloud  extinction  profile.  A  non-constant 
cloud  profile  will  have  an  effect  on  the  product  of 
cumulative  forward  scattering  and  depolarization 
ratio  and  will  most  likely  make  droplet  size 
recovery  unreliable. 


REFERENCES 

[1]  G.  Roy,  L.  R.  Bissonnette,  C.  Bastille,  and  G. 
Vallee,  1999:  Retrieval  of  droplet-size  density 
distribution  from  multiple  field-of-view  cross- 
polarized  lidar  signals,  Appl.  Opt.  38,  pp.5202- 
5211. 

[2]  I.  Veselovskii,  M.  Korenskii,  V.  Griaznov,  D. 
N.  Whiteman,  M.  McGill,  G.  Roy,  L.  Bissonnette 
“Information  contain  of  data  measured  with 
multiple-field-of-view  lidar,”  Appl.  Opt.  45,  6839- 
6848(2006) 

[3]  K.  Dukhyeon  and  J.  Lee,  2014:  Measuring 
cloud  droplet  effective  radius  and  liquid  water 
content  using  changes  in  degree  of  linear 
polarization  along  cloud  depth,  Opt.  Lett.,  39,  pp. 
3378-3381. 

[4]  D.R  Donovan,  H.  K.  Baltink,  ,  J.S.  Henzing, 
S.R.  de  Roode  and  A.P.  Siebesma,  2015,:  A 
depolarization  lidar-based  method  for  the 
determination  of  liquid-cloud  microphysical 
properties,  Atmos.  Meas.  Tech..  8,  pp. 237-266. 

[5]  K.  Sassen,  H.  Zhao,  1995:  Lidar  multiple 
scattering  in  water  droplet  clouds:  Toward  an 
improved  treatment,  Opt.  Rev.,  2,  pp.  394-400. 

[6]  G.  Tremblay,  G.  Roy,  2015:  High  fidelity 
imaging  algorithm  for  the  Undique  Monte  Carlo 
simulator.  27rd  International  Laser  Radar 
Conference,  New  York,  USA. 


